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We analyze the dynamics of a two-level system subject to driving by large-amplitude external 
fields, focusing on the resonance properties in the case of driving around the region of avoided level 
crossing. In particular, we consider three main questions that characterize resonance dynamics: 
(1) the resonance condition, (2) the frequency of the resulting oscillations on resonance and (3) 
the width of the resonance. We identify the regions of validity of different approximations. In a 
large region of the parameter space, we use a geometric picture in order to obtain both a simple 
understanding of the dynamics and quantitative results. The geometric approach is obtained by 
dividing the evolution into discrete time steps, with each time step described by either a phase 
shift on the basis states or a coherent mixing process corresponding to a Landau-Zener crossing. 
We compare the results of the geometric picture with those of a rotating-wave approximation. We 
also comment briefly on the prospects of employing strong driving as a useful tool to manipulate 
two-level systems. 



I. INTRODUCTION 

Two-level systems are ubiquitous in various fields of 
physics. A large number of quantum phenomena rely on 
the existence of two quantum states, or their underly- 
ing principles can be understood using the simple model 
of a two-level system. Recently, two-level systems have 
gained renewed attention as they represent the building 
blocks for quantum information processing (QIP) appli- 
cations 

In the study of two-level systems, as well as many other 
quantum systems, avoided level crossings are associated 
with a wide variety of interesting phenomena. Large 
amounts of literature have been devoted to analyzing the 
dynamics of a two-level system driven around an avoided 
crossing, particularly in connection with Landau-Zener 
(LZ) physics 0. These avoided crossing regions have a 
special significance in QIP applications because coher- 
ence times are usually longest in those regions, hence the 
term optimal point. Needless to say, nontrivial evolution 
of quantum systems is usually associated with some kind 
of energy-level crossing. 

In this paper we discuss the dynamics of a two-level 
system driven by strong ac fields in the vicinity of an 
avoided crossing. As mentioned above, numerous studies 
have been devoted to this problem, approaching it from 
different angles and applying it to different physical sys- 
tems. Here we present a geometric picture that leads to 
a simple understanding of the behaviour of the system. 
This approach can also be used to derive quantitative re- 
sults in regions where other methods fail. We also keep 
in mind the idea of trying to find useful applications for 
strong driving as a tool to manipulate two- level systems, 
or as they are referred to in the context of QIP, qubits 
(in this context, see e.g. 0,13, 



A good example displaying the richness of strongly 
driven two-level systems is the so-called coherent destruc- 
tion of tunnelling (CDT) Q- A particle in a symmet- 
ric double- well potential will generally oscillate back and 
forth between the two wells. If we now add an oscillating 
energy difference between the two wells, the frequency 
of the tunnelling oscillations changes. At certain combi- 
nations of the driving parameters, the tunnelling oscil- 
lations are frozen. This phenomenon has been analyzed 
from several different perspectives d, H, EL EH EH • 
We shall show below that it can be understood rather 
easily using our approach. 

Here we consider the situation where we start with an 
undriven system at an arbitrary bias point. We then an- 
alyze the dynamics that results from the application of 
strong driving fields. As one would expect, resonance 
peaks occur at properly chosen values of the parame- 
ters. We analyze the resonance conditions using two 
approaches: one employing a rotating-wave approxima- 
tion (RWA) and one employing an approximation of dis- 
cretized evolution characterized by a sequence of fast LZ 
crossings. Between the two approaches, and the well- 
known weak-driving case, most of the parameter space is 
covered. 

This paper is organized as follows. In Sec. II we present 
the model system, the Hamiltonian that describes it and 
some preliminary arguments. In Sec. Ill we present a 
RWA that can be used to describe the dynamics in a 
certain region of the parameter space. In Sec. IV we 
present a geometric picture that is useful to describe the 
dynamics in another region of the parameter space (note 
that there is some overlap between the validity regions of 
Sees. Ill and IV). Section V contains a discussion of the 
results and some concluding remarks. 
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II. MODEL SYSTEM AND HAMILTONIAN 

We consider a two-level system described by the Hamil- 
tonian: 



Hit) = - — a z- 



(1) 



where A is the (time-independent) coupling strength be- 
tween the two basis states [l4j; e(t) is the time-dependent 
bias point; and a x and a z are the usual x and z Pauli ma- 
trices, respectively Note that we take h — 1 throughout 
this paper. For definiteness and simplicity in the algebra, 
we assume harmonic driving, i.e. we assume that e(t) can 
be expressed as 



e(i) = e + A cos(ujt + . 



(2) 



where eo is the dc component of the bias point; and A, 
w and <p are, respectively, the amplitude, frequency and 
phase of the driving field. We shall, with no loss of gen- 
erality, take all the parameters in Eqs. fl} and ||3J| to be 
positive. In order to simplify the appearance of the ex- 
pressions below, we shall take = 0. This only simplifies 
the intermediate steps of the algebra, but it does not af- 
fect any of the main results. The energy- level diagram 
and the applied driving field are depicted in Fig. 1. 
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FIG. 1: (color online) (a) Energy-level diagram E(e), in blue, 
of a two level system with minimum gap A and (b) the time- 
periodic bias e(t), in red. The vertical dashed line represents 
dc bias point eo- 

We should mention here that recently this setup was 
studied experimentally using superconducting qubits in 
Refs. H, Uj, both of which contain theoretical analysis 
that overlaps with ours. Related theoretical studies can 
also be found in Refs. 0, Q. 

In the absence of driving, the behaviour of the system 
is simple. In the unbiased case (i.e., when eo = 0), the 
eigenstates of the Hamiltonian are the symmetric and 
antisymmetric superpositions of the states ||) and 
As a result, if the system is initially in one of these two 
states, it oscillates back and forth between them. In the 
strongly biased case (i.e., when e >• A), the states ||) 
and ||) are, to a good approximation, the eigenstates of 



the Hamiltonian. If the system is initially in one of them, 
it will only experience small oscillations, occupying the 
other state with a maximum probability (A/eo) 2 . One 
could say that in this case the oscillations, which are 
driven by the a x term in the Hamiltonian, are suppressed 
by the energy mismatch eo between the states ||) and ||). 

In the weak driving limit (i.e., when A <C -\/A 2 + ef), 
Rabi-oscillation physics applies. Resonance occurs when 
and the frequency of Rabi oscillations is 



= 



given by f2 = (Asina)/2, where the angle a is defined by 
the condition tana = A/eo- Rabi oscillations occur be- 
tween the eigenstates of the Hamiltonian (excluding the 
driving term), such that in the unbiased case the oscil- 
lations occur between the symmetric and antisymmetric 
superpositions of the states ||) and |J.). Higher-order 
processes, with nui — \/ A 2 + eg, can be described easily 
as well in this limit. We shall collectively refer to the 
weak-driving limit using the easily recognizable name of 
Rabi-physics limit. 

Since Rabi physics is well known [l5| , we shall not dis- 
cuss it in any detail here. Instead, we focus on the case 
w here the amplitude A is comparable to or larger than 
VA 2 + e 2 . A suitable RWA can be used to obtain a 
good description of the dynamics when lo 3> A. This ap- 
proach, described in Sec. Ill, can explain some interesting 
features of the problem, in particular the phenomenon of 
CDT mentioned in Sec. I. 

A different approximation can be used when (^4— eo) 3> 
A and Auj ^S> A 2 (note that this condition overlaps with 
the validity condition of the RWA) . In this case we can 
think of the dynamics as being composed of a sequence of 
LZ crossings separated by periods of free evolution of the 
basis states. We shall take this approach to analyze the 
problem in Sec IV. We shall then compare the validity 
conditions of the different approximations in Sec. V. 



III. HIGH-FREQUENCY DRIVING: 
ROTATING- WAVE APPROXIMATION 

We now take the system described by the Hamiltonian 
in Eq. (TTJ) and make a transformation to a rotating frame, 
such that a wave function 1-0) in the lab frame can be 
expressed as 



(3) 



where 



U(t) =exp|^ fe t + ^sinut\ ct z j , (4) 

and |^') is the wave function in the rotating frame. In 
other words, we use the interaction picture with the cou- 
pling term treated as a perturbation. The Schrodinger 
equation 



i-\i))=H{tm 



(5) 
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can now be written as 



We now make use of the relation 



with 



i^W) = H'{t)W), 



(6) 



H'{t) = U\t)H{t)U{t)-iU\t) 



dU(t) 
dt 







exp{izsin7} = Jn(z) e m7 , 



(8) 



n— — oo 



. v I .(7) find that 



where J n (x) are Bessel functions of the first kind, and we 



J 



E 4 



.4 



-i(nct;+eo)t 



(9) 



r 



Note that in going from Eq. JTJ) to Eq. ^ we have not 
made any approximations. It will be useful below to use 
the following asymptotic behaviour of the Bessel func- 
tions: 



Jn{z) 
Jn(z) 



2 n nV 
~2~ 



cos 



z- (2n + l)- 



z < 1 



z > n. (10) 



We now perform a rotating wave approximation 
(RWA): we assume that all the terms in the sum in 
Eq. @ oscillate very fast compared to the timescale of 
the (coarse-grained, or smoothened) system dynamics, 
with the exception of one term that is identified as the 
resonant or near-resonant term. The effect of the non- 
resonant terms can therefore be neglected when studying 
the long-time system dynamics. We should note here 
that this RWA is different from the one widely known 
and used in the case of weak driving. In that case the 
linear driving field is replaced by a rotating field, clock- 
wise or counter-clockwise depending on conventions. The 
other component of the field, i.e. the one rotating in the 
opposite sense, is neglected [171 ]. 

The resonance condition is identified using the intu- 
itive requirement that the resonant term of the RWA is 
constant in time. We therefore have the resonance con- 
dition: 

nu + e o = (11) 

for some integer n. With parameters satisfying the res- 
onance condition, the ||) <-> |J.) oscillation frequency is 
given by: 



fi = A 



I I- 

ui 



(12) 



One usually identifies the resonance with a given value 
of n as describing an |n|-photon process. The Rabi res- 
onance condition corresponds to the case n = — 1, where 
we find that u = eo (the difference from the condition 
uj = y/ A 2 + £q will become clear shortly). Assuming 
weak driving, we find that the frequency of Rabi oscilla- 
tions is given by: 



o A A 
Q = — x — . 

e 2 



(13) 



For large values of eo, the factor A/eo is a good approxi- 
mation to the factor sin a, where a is the angle between 
the static bias field and the driving field. The Rabi fre- 
quency vanishes asymptotically as eo — > oo, as it should. 
Note also that fl (Eq. fTS")) increases with increasing A for 
small values of A, but generally decreases as 1/ \[~A for 
large values of A. The mechanism responsible for this lat- 
ter behaviour will become clear in Sec. IV. In particular, 
note that the turning point between the two behaviours 
(or, alternatively, the maximum in f2) occurs when A is 
comparable to eo- 

The width of the resonance can be obtained from 
the following considerations. If the driving frequency is 
shifted from exact resonance (Eq. ITTj) by 5lo, the reso- 
nant (i.e., slowest) term in Eq. ([9]) oscillates with fre- 
quency nStu. When these oscillations become faster than 
the |t) *-* ||) oscillation dynamics, which is character- 
ized by f2, the resonance is clearly lost. The width of the 
resonance can therefore be taken as: 



(14) 



Using higher-order processes (i.e., with |n| > 1) therefore 
results in resonances that are narrow compared to the 
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on-resonance oscillation frequency. This property can be 
useful, for example, in applications where one is dealing 
with several closely spaced resonances. If one is trying to 
drive only one of those resonances, this approach provides 
a possibility to target a single resonance, without neces- 
sarily making the oscillation dynamics extremely slow. 

We now note that for the above resonance condition 
(Eq.QTJ) to hold we must be able to neglect all the oscillat- 
ing terms. This means that we require that |nw+£o| ^ A 
for all n except the one satisfying the resonance condition 
[l8| . Keeping in mind that \nto + e | = for one value 
of n, we find that the validity condition can be expressed 
simply as w > A. This explains the difference between 
the above resonance condition and the usual resonance 
condition for Rabi oscillations; for example, when eo = 
and ijj = A we cannot keep only one term in Eq. §§§ , and 
this approximation breaks down. Note, however, that 
even if eo — we can still use this approximation, as long 
as w > A. This approach is very useful when discussing 
the phenomenon of CDT, which we do next. 

The much-analyzed phenomenon of CDT Q can be ex- 
plained using the above approach. In the unbiased case 
(i.e., when eo = 0), Eq. (fTTj) gives the interesting result 
that, regardless of the value of u>, taking n — always 
satisfies the resonance condition. This means that oscil- 
lations between the states ||) and ||) will always occur 
with full \]) <-* | J.) conversion, provided of course that 
we can neglect all the terms with n ^ (or in other 
words to 3> A). This statement is obvious for no driving 
(A = 0), but it is not obvious that for large driving helds 
(A ^ A) full conversion should occur. From this point of 
view, it looks more surprising that full oscillations occur 
at all for strong driving, even though the s yst em hardly 
spends any time in the degeneracy region [19[ . Accept- 
ing the existence of these oscillations, we now take the 
oscillation frequency as given by Eq. (fl~2|) with n = 0. It 
is now clear that CDT occurs when 



■/„ [ - 



0, 



(15) 



such that the resonant term in Eq. ([9]) has a vanishing co- 
efficient (note that Eq. (fTS"]) agrees with the results of pre- 
vious work 0, Q). With parameters satisfying the CDT 
condition, the frequency of |t) <-> ||) oscillations van- 
ishes, and the oscillations are consequently suppressed. 
One could therefore say that CDT is simply the state- 
ment that the ||) «-> ||) oscillation frequency becomes 
extremely small if we use driving parameters that are 
close to a point satisfying Eq. (|15p . 



IV. REPEATED TRAVERSALS OF THE 
CROSSING REGION: TRANSFER-MATRIX (TM) 
APPROACH 



using the RWA above is still valid, that approach becomes 
less intuitive. Instead, one can gain better insight into 
the problem by analyzing the dynamics as composed of 
discrete steps, as we shall do in this section (see also 
Ref. i). 

Let us take the limit where (A — eo) S> A. We can now 
think of the system as undergoing a sequence of LZ cross- 
ings (represented by the red and green lines in Fig. 2). At 
each crossing event, the states |f) and | J.) experience some 
mixing. If we approximate the sweep across the degen- 
eracy region by a linear ramp of the bias point between 
two points located symmetrically around the degeneracy 
point (denoted by t c — r and t c +T, where t c is the time at 
which the bias point is at the center of the avoided cross- 
ing), we find that the LZ crossing can be approximately 
described by the evolution matrix: 



G 



LZ,k 



cos ■ 

. X t 
sin — e 
2 



X 



■ X 
sin — e 
2 

cos ■ 



>0l 



(16) 



where the angle \ is defined by the LZ transition proba- 
bility 



. aX -i / ^A 2 

sin — = 1 — exp < — - — 
2 v I 2v 



(17) 



v is the sweep rate, and k defines the direction of the bias 
sweep across the degeneracy region: k = 1 when e(t) goes 
from positive values to negative values and k = 2 when 
e(t) goes from negative values to positive values (Note 
that the sweep rates in the two directions are equal, and 
thus x is independent of k). The angles 6*lz k are given 
by 



0LZ,1 = 7T — 6*s 



tokcs 



l - f ^A 2 + (A 2 - e 2 )uj 2 T 2 dT 



7T - Stokes - 2^V A2 + - e 0^ V2 



A 2 



sinh 



7T - 6>stokcs - -\l A 2 ~ elujT" - 



A 2 



?LZ,2 



1 

" 2 

A 2 



log 



2VA 2 



s..,k,., - -\i A 2 - e 2 Q ujT 2 + 



2 V /I 2 ^4 i , 



log 



2y/A 2 - e 2 LOT 

A 



Wcs = -+arg[r(l-W)]+<5(ln£-l) 
5 = 4JT' 



(18) 



We now focus on the case of strong driving where the 
system repeatedly traverses the crossing region [i.e., (A — 
eo) ^> A]. Even in the case u> ^> A, where the treatment 



and T(x) is the gamma function (see e.g. Ref. [2(| for a 
detailed derivation). Note that the Stokes phase 6>stokcs 
approaches zero in the slow-crossing limit (i.e., when 
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6 — ► oo), and it approaches 7r/4 in the fast-crossing limit 
(i.e., when 5 — >■ 0) 22]. The description of LZ crossing 
processes using unitary matrices of the form of Eq. (fl"6|) 
is sometimes referred to as the transfer-matrix method, 
which wc follow here (the term scattering matrix can also 
be found in the literature). 

Between each two LZ crossings, the system moves far 
from the degeneracy region and acquires a relative-phase 
factor between the states |f) and ||), but with no mixing 
between them. Because of the asymmetry (i.e., the fact 
that, in general, e ^ 0), there are two phase factors 
corresponding to the system being on the right or left 
side of the degeneracy region. We therefore find that 
between crossings, the system evolves by the evolution 
matrices: 

4 = ( e 7 X ) * ^ 

where j represents the two sides of the degeneracy point 
(we shall refer to them as 1 and 2). 

Before starting to analyze the dynamics resulting from 
combining the above matrices, we note that the values 
of the angles 6\, 82, 0lz,i and #lz,2 depend on where 
we set the boundaries between the different time steps 
(i.e. they depend on r). It is therefore useful to define 
the boundary-independent phase factors Q\, 62, &lz,i and 
0lz,i- We define 6\ and 62 as half the relative phases ac- 
cumulated between the energy eigenstates between the 
times of two successive level crossings (ignoring in this 
evaluation the fact that there is mixing dynamics occur- 
ring between the states): 

§! = - J ^V( e o + A cos cot) 2 + A 2 dt 

= - J ^(eo + Acosujt)dt- f x (e ,A,uj, A) 

VA 2 - 4 eo _i -eo , , „ a \ 
= 1 cos — — - /i(e , A,uj, A) 

LO LO A 

VA 2 - eg e _! e 7re , , . . > 
= - + — cos L — fx{e ,A,u,A), 

OJ U! A LO 

- 2 = V^ 2 - 4 _ £0 cqs -x eo + M ^ A) (20) 

LO LO A 

where f\ and f% are functions that, roughly speaking, are 



proportional to A 2 /(Alo) times the logarithm of A/ A (see 
Eq. [TH]). With this definition of 0\ and 62, we find that 
the appropriate values of #lz.i and 0lz,2 are given by: 

#LZ,1 = 7T — ^stokes 

9hZ,2 = ^Stokes- (21) 

We now have a sequence of finite time steps, each of 
which is described by a simple evolution matrix, as shown 
in Fig. 2. The dynamics over a large number of driving 
cycles can be constructed by multiplying the evolution 
matrices. The evolution of the system over N cycles is 
e(t) I 




FIG. 2: (color online) The time evolution of the system di- 
vided into discrete time steps, each of which can be described 
by a simple evolution matrix. 



therefore described by the matrix 

Gn — ( Glz,2G2Glz,iGi J , (22) 

with possibly some minor modifications at the beginning 
and end of the sequence, depending on the exact initial 
and final bias points [2l[ . We now take a single cycle and 
calculate the evolution matrix that describes it. We find 
that 

Glz,2 G2 Glz,i G\ = ( ^l 2 ) , (23) 

V 9i2 9 11 J 

where 



gn = cos 2 I e~^ 1+e2) - sin 2 | S^, 1+ e^e 1+ e 2 ) 

gi2 = sin I cos I ( e i(«Lz,i+9i-fa) + e r(e^+e 1+ e^ (24) 
It is useful to rewrite Eq. ([23]) as 

sin^T e^ FC ^ / e -^c/2 \ 
G LZ ,2 G 2 G LZ ,i Gi = j r _ 2 J ^ I y Q e , eF U c /2 J , (25) 
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where the subscript FC indicates that the above equation describes the evolution of the system over a full cycle of 
the driving field. We now find that 

• 2 CFC . . 2 X 2 ( 0LZ,l - @LZ,2 

sin — — = 4 sin — cos ' 



2 2 V 2 

cos 2 x/2 sin(6»i + 6> 2 ) + sin 2 x/2 ma(9 LZ ,i - 6 LZ . 2 + 6 1 - 9 2 ) 



9pc = 2 arctan ■ 

cos' x/^cos^i + v 2 ) - sin" x/' ZC0S {Vlz,i - Vlz.i + Vi - V2) 

, 0LZ,1 + QhZ.2 + 2^1 — 6>FC ,„,,>. 

he = ^ ■ (26) 



We can now analyze the dynamics described by 
Eq. ([25]) . The first matrix in the product describes a 
rotation by a some angle Cfc around an axis in the xy- 
plane. The second matrix in the product describes a 
rotation by an angle $fc around the z axis. Since we 
are looking for resonance-like dynamics (which naturally 
implies the existence of a resonance condition), we want 
the angle £fc to be small. This condition is satisfied in 
two opposite limits: the fast- and slow-crossing limits. 
We shall treat these two limits separately below, after 
making the following general observations. 

With the above geometrical interpretation of the roles 
of the angles £fc and Afc, the resonance condition is 
clear. If #pc is a multiple of 2tt, the z-axis rotation does 
not affect the dynamics, and the small rotations of (fc 
add up to produce full oscillations between the states ||) 
and ||). If, on the other hand, 6*fc takes a value that 
is different from any multiple of 2tt by more than £fCj 
the small rotations of £pc wm n °t add up in an ideal 
manner, and the oscillations will be suppressed [23[. In 
the following subsections we treat the two limits where 
simple analytic formulae can be obtained. 

The width of the resonance can also be obtained using 
the geometrical picture explained above. We require that 
the phase factor Ope be within a distance £fc from a 
multiple of 2ir. Starting from Eq. ([50]) . and writing 



#FC 



2lTe ° A 
Zirn — ouj, 



(27) 



where Suj is the deviation from exact resonance, we find 
that the width of the resonance is given by: 



5u) 



^ 2 Cfc 
2vre 

n 



(28) 



As in Sec. Ill, we find that in addition to the usual fac- 
tor of oscillation frequency, the width of the resonance 
now contains the factor 1/n. This means that with the 
proper choice of parameters, the width of the resonance 
can be made substantially smaller than the on-resonance 
oscillation frequency. 



A. Fast-crossing limit 

In this subsection we assume that each crossing is tra- 
versed in the fast limit. This means that we require the 
sweep rate across the degeneracy region (i.e., lj\J A 2 — e 2 ,) 
to be much larger than the square of the gap size of the 
crossing A 2 . The LZ probability in this limit is given by 



2 X 



7rA 2 

~2v 



(29) 



Note that in this limit sin 2 (x/2) < cos 2 (x/2). 

The resonance condition for the constructive accumu- 
lation of small rotations can now be obtained by noting 
that: 



e FC w 2(e 1 + e 2 ) 

rT+2-rr/uj 



dt [eo + A cos ujt] 



27re 



The resonance condition is therefore given by 



— = n 

UJ 



(30) 



(31) 



for some integer n. This is the same condition that we 
found in Sec. Ill, using an entirely different approach 
(note also that the validity conditions of the two ap- 
proaches are different, a point to which we shall come 
back in Sec. V). The angle </>fc is given by 



'FC 



7T 

2 



(32) 



Using the approach of this section it might seem some- 
what surprising that the amplitude A does not appear 
in the resonance condition, even though A can be the 
largest energy scale in the problem. One should also 
note here that Eq. (fTTj) is only approximate in the weak- 
driving limit, with a correction that is proportional to A 2 
(i.e. the Bloch-Siegert shift). One might therefore expect 
that such corrections will take over at some point, such 
that the amplitude A becomes an essential part of the 
resonance condition. Counterintuitively, however, these 
corrections change behaviour and vanish asymptotically 
for large A (in the fast-crossing limit). The fact that the 
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two (identical) resonance conditions are good approxima- 
tions deep in opposite limits demonstrates further that 
the derivation of Eq. (|3"T|) in this section should not be 
thought of as simply a re-derivation of Eq. (fTTj) . It is 
worth noting here that even for non-harmonic driving we 
can follow a similar analysis to what was done in this 
section and find that the resonance condition is still in- 
dependent of the driving amplitude, with eo replaced by 
the time-averaged value of the bias point and u> replaced 
by 2tt over the driving period [see Eq. ([3TJ]) ]. 

The frequency of oscillations on resonance can now be 
obtained rather straightforwardly. If we assume that the 
resonance condition in Eq. (|3ip is satisfied, we find that 
the |t) <-> ||) oscillation frequency is given by: 



n = 



u>Cfc 
2tt 

2^ . X 
— sin — 

7T 2 



We therefore find that 



'LZ,1 



?LZ,2 



n 



2lu 



ttA 2 



7T V 2uy/A 2 -(% 



cos ( 9 



(33) 



(34) 



In the fast-crossing limit, the functions fx and f% in 
Eq. (|2"D|) are negligibly small, and the phase 6 2 is ap- 
proximately given by 



\A4 2 - 4 £ o _i eo 
— cos — . 

UJ U! A 



(35) 



In the special case eo — 0, we find that 62 — A/u>, and we 
recover Eq. ([P2]) with n = for the oscillation frequency. 

One might wonder why several expressions above are 
not symmetric with respect to 9\ and 9 2 - This asymme- 
try results from our grouping of evolution matrices into 
full cycles [see Eq. (|22pj. as well as the order of matri- 
ces in Eq. (|25p . These are clearly matters of convention. 
Different orderings of the matrices can result in different- 
looking expressions. The end result must of course be 
independent of this choice. For example, if we substi- 
tute Eq. (|33|) into Eq. ([3"4"]) , we do not see any convention 
dependence. 

In order to illustrate the transfer-matrix picture, we 
show in Figs. 3-5 numerical simulations of the dynam- 
ics in the validity region of that picture. We plot the 
occupation probability of the state |f) as a function of 
time, assuming that the system was initially in the state 
We can see in Fig. 3 that the occupation probabil- 
ity exhibits sudden jumps that correspond to LZ cross- 
ings. The steps are rather large in this figure because 
the crossings are not in the fast limit. If we look on long 
time scales, we can see that the dynamics looks like si- 
nusoidal oscillations. This long-time behaviour becomes 
particularly smooth when the transition probability in a 
single LZ crossing is small, as is the case in Fig. 4. We 
also plot in Figs. 4 and 5 sinusoidally oscillating functions 



with frequencies given by Eq. ([34]) [24|, as well as sinu- 
soidally oscillating functions with frequencies given by 
Eq. (fl"2"]) from Sec. III. In Fig. 5 the driving frequency u> 
is smaller than A. We therefore find inconsistency in the 
predictions of Eq. (|12p . As a general rule, the RWA gives 
reasonable or good agreement with the numerical sim- 
ulations when the resulting oscillation frequency in the 
system dynamics is large. When the oscillation frequency 
is small, the effect of the resonant term is not necessarily 
large compared to that of the other terms in Eq. and 
the RWA fails. This is most clearly seen in Fig. 5(c). In 
Fig. 4, we are in the region where uj > A, and we always 
find that both Eq. ([12]) and Eq. ([34]) agree well with the 
numerical simulations. It should also be noted that since 
Fig. 5 corresponds to parameters that are not deep in the 
fast-crossing limit, Eq. (|34p shows some deviation from 
the true oscillation frequency (see Fig. 5a). 



pT 0-5 




FIG. 3: (color online) Occupation probability Pf of the state 
j|) as a function of time (in dimensionless units), assuming 
that the system was initially in the state |j). We take to/ A = 
5, and uj/A — 1. The driving amplitude is given by A/ A — 30 
for the black line and A/ A = 34.95 for the green (gray) line 
[note that the resonance condition is satisfied in both cases]. 
The transfer-matrix method and the RWA of Sec. Ill both 
agree well with the numerical results (their predictions are 
not shown here for clarity). 



The fact that the resonance condition is always satis- 
fied in the unbiased case is very clear in this approach. 
In this case, the phase factors 6\ and 9 2 accumulated on 
the two sides cancel because of symmetry, regardless of 
the driving amplitude and frequency. 



The factor | cos($2 — 7r /4)| in Eq. ([M]) gives a non- 
trivial dependence of the ||) «-> ||) oscillation frequency 
on the bias and driving parameters. It lies behind the 
phenomenon that even if the resonance condition is sat- 
isfied, it is still possible for the oscillations to be so slow 
that the resonance is effectively destroyed (see Fig. 3). In 
other words, it describes the same mechanism responsible 
for CDT, and it agrees with the Bessel ladder structure 
discussed in Ref. [3|. 
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At 

FIG. 4: (color online) Occupation probability Pj of the state 
||) as a function of time (in dimensionless units), assuming 
that the system was initially in the state ||). In all the fig- 
ures eo/A = ui/A = 3. The driving amplitude is given by 
A/ A =10(a), 15(b) and 20(c). The blue dashed line gives the 
(coarse-grained) predicted dynamics from the transfer-matrix 
method in the large-amplitude limit, and the red dotted line 
gives the predicted dynamics from Sec. III. 




At 

FIG. 5: (color online) Occupation probability P-\ of the state 
j|) as a function of time (in dimensionless units), assuming 
that the system was initially in the state ||). In all the figures, 
eo/A = 1, and to/ A = 0.5. The driving amplitude is given by 
A/ A =12(a), 16(b) and 20(c). The blue dashed line gives the 
(coarse-grained) predicted dynamics from the transfer-matrix 
method in the large-amplitude limit, and the red dotted line 
gives the predicted dynamics from Sec. III. 



B. Slow-crossing limit 



When A 2 /uj v / A 2 - e 2 , » 1, we find that 

cos -«exp|- — 

and sin 2 (x/2) > cos 2 (x/2). 
The angle #fc is now given by 



-2tt- 
+2/i 



2ire A JW 



4 — cos 



2/2, 



1 fo 

A 



(37) 



7FC 



?LZ,1 — PLZ.2 i- Cl — 02) 



where fi and f 2 grow, roughly speaking, as 
( 36 ) (A 2 /Aj) log(A/A) (see Eqs. M and [20j). Because 
there are no simple expressions for /i and f% (combined 
with the fact that they are much larger than unity in the 
slow-crossing limit), deriving an analytic expression for 
the resonance condition is not as straightforward in this 
case as in the fast-crossing limit. The resonance lines can 
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be obtained numerically, since a numerical calculation 
of 9i and 62 would be straightforward. Alternatively, 
one can obtain a rough idea about the shapes of the 
resonance lines by ignoring /1 and ji in comparison to 
the other terms in Eq. (|37p . Equating 6>fc to 27m, the 
rough approximation of the resonance condition with 
the above simplification is given by 



£0 



iru> 



2 — cos — 

7ru; A 



n. 



(38) 



In the A-€q plane, the resonance lines form arcs around 
the origin (i.e. around the point A = e = 0). Note that 
when eo = 0, the above resonance condition reduces to 



2A 

TTUJ 



(39) 



such that the parameter A is crucial in the resonance 
condition, in contrast to the results of the fast-crossing 
limit. The oscillation frequency can also be calculated 
numerically using Eq. (|33p . The angle </>fc is given by 



DISCUSSION AND CONCLUSION 
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FIG. 6: (color online) Regions of validity for different approxi- 
mations. TM stands for transfer-matrix method, RWA stands 
for the rotating-wave approximation presented in this paper, 
and Rabi stands for the weak-driving limit best-known in con- 
nection with Rabi oscillations. The axes are the frequency 
uj and amplitude A of the driving field, both normalized to 
the minimum gap A. The Rabi region is described by the 
condition A/ A < 1 (regardless of the value of ui), and it is 
shown by green squares. The RWA region is described by 
the condition u/ A > 1, and it is shown by red ellipses. The 
TM region is described by the conditions A/ A > 1, and it is 
shown by blue symbols. The TM region can be divided into 
three sub-regions depending on the parameter Aui/A 2 : the 
slow-crossing limit (x symbols), the fast-crossing limit (cir- 
cles) and the intermediate-speed case (stars). Note that eo 
was generally assumed to be comparable to A in this figure. 



In this paper we have presented two approaches to 
study the problem of a strongly driven two-level sys- 
tem. The first one (presented in Sec. Ill) was based on 
a rotating-wave approximation, whereas the second one 
(presented in Sec. IV) was based on a discretized descrip- 
tion of the dynamics. It is important to note that the two 
approaches have different regions of validity, as shown in 
Fig. 6. We also include in Fig. 6 the region of validity 
for the weak-driving limit, which is most clearly associ- 
ated with Rabi-oscillation physics. As can also be seen 
in Fig. 6, the different regions overlap substantially, and 
two of them can sometimes be used to describe the same 
situation. For example, although the derivation and the 
appearance of the results of Sec. IV rely strongly on the 
picture of LZ crossings, the results agree with those of 
Sec. Ill whenever w> A. Naturally, the deeper one goes 
into one of these regions, the better the results one can 
expect to obtain from the corresponding approximation. 

The TM region is divided into three sub-regions. When 
Auj ^> A 2 , each LZ crossing occurs in the fast limit and 
the expressions given in Sec. IV. A are valid. In the op- 
posite limit, i.e. when Auj <C A 2 , the crossings occur in 
the slow limit, and Sec. IV. B applies. Between these two 
limits, one has the intermediate-speed, or general, case. 
Even though we have not derived any quantitative results 
describing the dynamics in this case, it can be conceptu- 
ally understood using the TM picture discussed here (as 
can be seen from Fig. 5). 

The approaches discussed here therefore cover a large 
portion of the parameter space. They provide alterna- 
tive points of view for understanding the mechanisms at 
play in the dynamics of this system (25j. It is worth 
noting that the approach of Sec. IV is not limited to har- 
monic driving. It can be used to treat any system with 
large-amplitude driving around the degeneracy point, as- 
suming the approximation of linear sweeps through the 
crossing region is valid [26| . 

Experiments on two-level systems have generally suf- 
fered from short coherence times. With the advent 
of the field of QIP, the need for long coherence times 
has spurred a fast advance in the direction of isolating 
qubits from their environments, thus resulting in rela- 
tively long coherence times. For example, high-order 
processes and quantum interference between LZ crossings 
have already been observed in superconducting qubit sys- 
tems p| El HI| ■ One could in the future realistically 
think about using strong driving as a tool to manipulate 
qubits. The mechanisms discussed in this paper can be 
used in constructing such qubit-manipulation tools. 
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